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Low temperature quantum corrections to the density of states (DOS) and the conductivity are 
examined for a two-dimensional(2D) square crystal with substitutional impurities. By summing the 
leading logarithmic corrections to the DOS its energy dependence near half-filling is obtained. It is 
shown that substitutional impurities do not suppress the van Hove singularity at the middle of the 
band, however they change its energy dependence strongly. 

Weak disorder due to substitutional impurities in the three-dimensional simple cubic lattice re- 
sults in a shallow dip in the center of the band. 

The calculation of quantum corrections to the conductivity of a 2D lattice shows that the well- 
known logarithmic localization correction exists for all band fillings. Furthermore the magnitude of 
the correction increases as half-filling is approached. The evaluation of the obtained analytical re- 
sults shows evidence for delocalized states in the center of the band of a 2D lattice with substitutional 
impurities. 
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I. INTRODUCTION 

The investigation of effects of weak localization in low-dimensional electronic systems still preserves its popularity 
today, since experimental studies of these systems give occasionally surp|rising results like metal-insulator phase 
transitions in the two-dimensional (2D)_electron gas at zero magnetic fieldtl and unusual behavior of dephasing time 
in quasi one-dimensional metallic wireso at low temperatures. Impurity effects in one-dimensional (ID) disordered 
systems have been studied accurately due to the existence of methods which give exact results for the cases of 
both weak and strong disorder .Uu Also perturbative approaches give good results for three-dimensional (3D) weakly 
disordered systems. 

Weak localization corrections to the kinetic coefficients of 2D disordered systems are logarithmically divergent 
corrections and it is hard to sum perturbatively the leading divergent contributions. The dimension two is the marginal 
dimension for the localization problem and a small external perturbation can change the character of localization in 
these systems. ^ 

According to the scaling theory of Abrahams et. ala all states of ID and 2D electronic gases moving in the field 
of randomly distributed impurities are completely localized irrespective of the degree of randomness, and in 3D the 
Anderson metal-insulator phase transition occurs with increasing impurity concentration. Notice that the result of 
the scaling theory for ID disordered systems is in agreement with exact results.Bfl The scaling theory, and also the 
diagrammatic approach to the problem for 2D weakly disordered systems has revealed logarithmic quantum corrections 
to the conductivity,!!! which tend towards localization. The logarithmic correction to the conductivity found in [ 
can be expressed as 

T* 

^-^o = -^ln- (1) 

where cto is the Drude conductivity, r* = min{T;^, ^} with r and being the elastic and inelastic scattering 
times, respectively; L is the linear size of a 2D system and to is the external frequency. 

Notice that weak disorder gives no singular contribution to the density of states (DOS) p(e) of electron gas models 
when the condition ppl ^ 1 or epT 3> 1 is satisfied {pp and ep are the Fermi momentum and the Fermi energy; I 
is the mean free path). AccordinE-ip the Einstein relation a = e^pD the quantum correction Eq.(^ to a is due to 
changes in the diffusion coefficientlDfl D. „ |— , 

Even short-range and weak correlations in weakly disordered metals have been shownQ~!12! to result in nontrivial 
quantum corrections to the DOS of low-dimensional conductors near the Fermi level. Electron-electron interactions 
in disordered systems also give rise to quantum corrections to the conductivity, which tend to localize the electronic 
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states. Further, the DDS_pf strongly doped semiconductors in the presence of long range Coulomb interaction vanishes 
at the Fermi surface,L50 which is commonly referred to as Coulomb gap. 

Disordered metals in most of the papers on localization problems are modeled as a free electron gas moving in 
the random field of impurities. However, doped crystals with low concentrations of impurities usually preserve their 
periodical structures and the impurity atoms in most cases substitute the host atoms of the lattice. Substitutional im- 
purities appear to have strong influence on the physical properties of low-dimensional lattices near the commensurate 
values of the electron wave length A with the lattice constant a [ |l^ . 

In this paper we shall study the effects of substitutional impurities with small concentration on the density of states 
and the conductivity of a 2D square lattice. The effects of commensurability on the DOS for the 3D simple cubic 
lattice will also be considered for the sake of completeness. 

Bragg reflections of the electronic wave on the Brillouin zone boundary in the process of impurity scattering become 
essential as the middle of the band is approached. This process introduces a new relaxation time for umklapp scattering 
to the problem, which increases for a deviation from half-filling due to the distortion of the commensurability condition 
for the Bragg reflection. 

The lattices under consideration have square structure for 2D and simple cubic structure for 3D crystals. Our 
calculations show that umklapp scattering strongly changes the energy dependence of the DOS close to half-filling, 
causing it to decrease. For the 2D square lattice the DOS of the pure system has a logarithmic van Hove singularity at 
the middle of the band. This peak is shown to be preserved for the lattice with substitutional impurities. However its 
energy dependence is changed strongly, and the peak becomes narrower .113 For the 3D lattice, the van Hove singularities 
lie far from the center of the band and umklapp scattering forms a shallow dip on the Fermi surface, the depth of 
which increases with the impurity concentration. 

Periodicity causes the appearance of a new class of diagrams which give a contribution to the conductivity cr{uj). 
These corrections to the conductivity can be separated into two classes: corrections due to the diffusion coefficient and 
due to the DOS. Although quantum corrections to the diffusion coefficient slightly increase cr(w), the total corrections 
sum up to a decrease of (t{u) due to DOS corrections in the vicinity of half-filling. 

The paper is structured in the following form. In Sec.^ the method is described, where new diffuson and cooperon 



blocks due to umklapp scattering are introduced and calculated. Sec. [II is devoted to the calculation of corrections 
to the DOS near half-filling for 2D and 3D lattices. In Sec.|v[ the calculation of the conductivity of the 2D square 
lattice with substitutional impurities is presented. In Sec.^, a conclusion and a discussion of the results are given. 

II. DESCRIPTION OF THE METHOD 

The Hamiltonian of a d-dimensional simple cubic lattice with substitutional impurities is given by 

H = Ho + V{r) (2) 

where Hq is the tight -binding Hamiltonian of noninteracting electrons in a regular lattice with lattice constant a and 
nearest-neighbor hopping; ^(r) = J^i U{r — Hi) is the impurity potential with R.^ being the positional vector of an 
impurity randomly located on the i-th lattice site. 

By introducing the second quantization operators ^ and Cp^a for an electron with momentum p and spin cr, 
Eq.(||) can be rewritten as 

H = '^e{p)cl,,Cp,„+ pi,„p(q)C/(q)cj,_„Cp+q+G,<T (3) 

where /3imp(q) = L~''-'^-exp{iqR,i) and U{q) is the Fourier transform of a single impurity potential. Throughout 
this paper, we work in units of ?i = 1. The momenta p and q vary in the first Brillouin zone and G is a reciprocal 
lattice vector. e(p) in Eq.(^) is the energy spectrum of an electron on a c?-dimensional square lattice, 

e(p) = t[2 — cos{pxa) — cos{pya)] for d = 2 (4) 

e(p) = t[3 — cos{pxa) — cos{pya) — cos{pza)] for d = 3 (5) 

where Px.y.z — ^]^^ — ^x,y,z with < Tir^ y ^ < ^''^'^ and t is the tunneling integral for nearest-neighbor sites. 

The electronic bandwidth is W — At for d = 2 and W — 6t for 3D systems. As half-filling is approached, which 
corresponds to ep — 2t for 2D and = 3i for 3D system, the Fermi surface becomes nested, i.e. it exhibits surface 
elements which can be mapped onto each other through one nesting vector. In our case, the nesting condition is 
satisfied for all parts of the surface, and for 2D square lattices, the Fermi surface at half-filling is flat (see Fig.|l|). 
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The second part of the Hamiltonian H in Eq.(||), which describes an electron scattering on randomly distributed 
impurities, contains both normal (for G = 0) and umklapp (for G Tt— D) processes. The impurity concentration 
is assumed to be small, so that the Born approximation is appropriatetj to estimate the scattering process on the 
impurities. The impurity potential in this case is chosen to be a (5-correlated Gaussian potential with zero average 

...... 

Conventional diagrammatic techniqueaiS are applied to calculate the effects of umklapp scattering on the DOS and 

the conductivity. The bare Green's functions G'5j^(e,p) at zero temperature are given as 

G°«,A(e,P)- , , , » • (6) 

e- (e(p) - ef) + 27Signe 

where e > and e < correspond to retarded G^j and advanced G\ Green's functions, respectively. The energy 
spectrum e(p) in Eq.(^ is expressed by Eqs.(^) and (||) for 2D and 3D cases, respectively. For momenta lying close 
to the Fermi surface (in the first Brillouin zone) the energy spectrum can be linearized around the Fermi surface and 
e(p) — ep K. Vi?(p — Pi?) = WFdpl — Vf) cos a, where v(p) = ^ = ai{sinp^a, sinp^^a} is the group velocity of the 
electron wave packet. In contrast to the electron gas model the linearized energy spectrum contains the additional 
parameter cos a due to non-coUinearity of the momentum and velocity vectors. For small band filling, the Fermi 
surface is very similar to a sphere, so this factor can be set equal to unity. Approaching half-filling, this parameter 
becomes weakly varied, e.g. < cos a < 1 for a 2D system at half-filling, where the range of variation of cos a 
is maximal. As an approximation, we can always set cos a = 1 . In the vicinity of commensurate points, especially 
at half-filling, electron scattering on impurities with large momentum transfer involves umklapp processes. The 
procedure of expansion of the energy spectrum around the Fermi surface can in this case be performed after the 
separation of the large momentum transfer Qq. 

In weak localization theory the maximally crossed diagrams are responsible for the low temperature quantum 
corrections to the kinetic properties of low-dimensional disordered systems, modeled as an electron gas moving in the 
field of randomly distributed impuritiesQ. These diagrams can be redrawn as ladder diagrams in a particle-particle 
channel (Fig.|a). The cooperon block has a pole for small total momentum of particles, |q|/ <C 1, and for small energy 
difference |w|t <C 1. The cooperon block turns out to have a diffusion pole also for large momentum transfer when 
the total momenta of the particles are (Qo + q) with |q|Z ^ 1 and when the total energy is small. In this case each 
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act of scattering on an impurities involves a large momentum transfer close to Qo- This process of course takes place 
mainly around half-filling when Qo = 'ipp, and it consists of simultaneous Bragg reflection of the electron in the 
process of scattering on the impurity, with relaxation time t-^. This scattering process will be further referred to as 
TT-scattering. 
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FIG. 2. The ladder series for a) the normal cooperon Cjv, b) the 7r-cooperon Ctt and c) the 7r-diffuson D^. In b) and c), Cr 
and D.,^ are expressed by — J^C^ and -D^ — J-^Dt^, where J.^ is given in Eq.dltl). 



Two different scattering processes and corresponding relaxation times can be separated: Normal scattering with 
T^v (see Fig.^), when the scattering on an impurity maintains the electron's momentum inside the first Brillouin 
zone after scattering and 7r-scattering, with a relaxation time (Fig.^). Here, momentum conservation is violated 
and this scattering process corresponds to an umklapp process. The total relaxation time r can then be expressed as 
i = + Notice here that the expression for the normal relaxation time tn according to Fig.|^ for the 2D case 
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can be given as 



- = WyJ ^|c^(pf,s)P (7) 

TN {2ttY Js |Vk| 

We assume tn to be constant and independent of the band filling. This means that the singularities of the integrand 
due to vanishing velocity at the saddle points are compensated by appropriate zeros of the impurity potential. 

To illustrate diagrammatically the impurity scattering including umklapp processes we represent the Green's func- 
tion of an electron with large momentum by a dashed line. The ladder series for the cooperon block with large 
momentum transfer, C7r(q, e), is diagrammatically shown in Fig|2|3. Notice here that the vertices c) and d) in Fig.^ 
are irrelevant for the cooperon blocks C7r(q, e) and CAr(q, e) (Fig. §a,b). Indeed, the dashing of the line alone has 
no physical meaning, it has only meaning for an electron scattering with large momentum transfer. In ladders, the 
vertices c) and d) of Fig.^ can always be transformed into normal vertices or those of Fig.^jb by proper redefinition of 
the internal momenta. 
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FIG. 3. Impurity vertices which give a contribution to a) the normal relaxation time tn, b) the relaxation time for 
TT-scattering r^. The vertices c) and d) do not give a contribution to vr-scattering. 

In the process of propagation in the particle-particle channel the momenta of the particles for the normal cooperon 
block Cjv(q, uj) lie on opposite sites of the Fermi surface during the diffusion process. In contrast to this, the momenta 
for the two particle propagator C-^ (q, e) with umklapp process lie on the same section of the Fermi surface and each 
act of impurity scattering coherently relocates the particles to the opposite section of the Fermi surface due to Bragg 
reflection. We will refer to this two-particle propagator Cr (q, e) as 7r-cooperon. 

The same situation occurs for the 7r-diffusion block DTr(k, e, oj) (Fig||) with a large difference (Qo + k) of electron 
and hole momenta, with |k|^ ^ 1. The comparison of the normal diffuson block D^Ck, e,uj) with the 7r-diffuson 
D^(k,e,Ci;) reveals again strong differences between these two processes. In the process of propagation in the TV- 
diffusion channel , electron and hole lie close to each other on the Fermi surface, with momenta p and p + k, and 
|k|Z ^ 1. They diffuse on the Fermi surface and keep the proximity of the momenta after each act of scattering. In 
contrast to this, the momenta of electron and hole for the 7r-diffuson block are located on opposite sides of the 
Fermi surface. Each scattering on the irnpurities interchanges the positions of electron and hole on the Fermi surface. 
Now let us sum the series shown in Fig^ for Cjv, Ct^, and Dt^. 

Summing up the ladder series in Fig.m, the following expression for Civ(q, t^^) is obtained: 

CAr(q,c.)^a^p|L/(0)p{ ,y^\~,'^!^'^'^if ' , e{-eie + u;)) + e{e{e + u;))} (8) 
^y/{l-iT\uj\)^ + {qiy - J 

Far from half-filled electronic band, ^ tn and r w tn. Therefore, CAr(q, has a diffusion pole for \llj\t ^ 1 
and \q\l <C 1. As half-filling is approached, umklapp scattering is intensified and ^ tn and t « t^. As a result. 
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the diffusion pole of C7v(q, e) disappears at half-filling. In this regime, the Fermi surface of the square or simple 
cubic crystal becomes nested (see Fig.|l| for the 2D case) with a nesting vector Qo = {±-^,±^} for the 2D lattice and 
Qo = {±-^,±-^,±^} for the 3D lattice. In this case, the following particle-hole symmetry of the electron dispersion 
with respect to the vector Qo holds for the half-filled band: 

e(P + Qo)-eF = -[e(p)-eF] (9) 

The main contribution to the low temperature properties now gives the 7r-cooperon C,r{q, e,Lu) in Fig.^. One gets 
the following expression for C\{q, e, us) by summing the ladder series in Fig.0: 



a(q,e,c.) = amp|C/|'Y^ (10) 



where Jtt is the expression for an elementary "bubble" in the ladder series and 



J.^Cin.pJ ^|C/pG(p,e),G(-p + q + Qo,e + c^) (11) 

To calculate Jtt, the large momentum Qo is removed using the electron-hole symmetry relation Eq.(^, and after this 
the energy spectrum is linearized around the Fermi surface. As a result we get the following expression for C^(q, e, lu): 

CM. = a..pl^(2p.)P{ _ ^^1^^ ^ _ 9(.(. + .)) + e(-.(. + .))} (12) 

So the TT-cooperon has a diffusion pole for total momenta close to Qo (|k|Z ^ 1) and small total energy |2e+a;|T ^ 1 
of the particles. The diffusion block _D^(k, e, w), which differs from C^(q, e,^) through time reversal of one electron 
line, has a pole for large (Qo + k) momenta difference with <C 1 and small total energy of electron and hole (see 
Fig.^). The calculation of Z?7r(k, e,w) is similar to that for C7r(q, e, t^-") and we obtain 

DAK^.^)-C,..,\U{2p,)\^{^-^^^ (13) 

After this, we can calculate the corrections to the DOS and the conductivity due to umklapp processes. 



III. DENSITY OF STATES AT HALF-FILLING 

The one-particle DOS of the regular d-dimensional lattice is expressed as 



where dS is an element of an isoenergetic surface in d-dimensional space, pl^ has a van Hove singularity at the points 
where the group velocity of the electron wave packet = Ve(k) vanishes.EZl 

The DOS of a clean ID lattice with constant spacing a is easily calculated to be pl^'' = ^ This expression 

Tza^ e{2t~e) 

shows that p^^^ is a regular function of the energy near half-filling when e ep — t and it has a power-like singularity 
when e approaches the band edges e ^ 0, 2t. The bare DOS of a 2D crystal with an energy spectrum given by Eq.(|^) 
has a logarithmic van Hove singularity in the middle of the band: 



dS 



Here and below, the electron energy e is measured from the middle of the band, i.e. e = corresponds to half-filling. 
The van Hove singularity of simple 3D crystals with nearest-neighbor hopping exhibits cusps near |e| — t: 

p^i^ = const - ^^^3^/2^3 ^^ (16) 
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In this section we will study effects of substitutional impurities on the DOS of 2D and 3D crystals. Randomly 
distributed impurities in d-dimensional electron gases have no effect on the DOS. 

The DOS of 2D and 3D simple (cubic) crystals with substitutional impurities turns out to have a quantum correction 
near the middle of the band even for noninteracting electrons. Effects of commensurability on the PjOS, and on the 
kinetics of ID disordered crystals near the middle of the band have been calculated by many authors. EJca Dyson first 
pointed outE£l that the DOS of phonons of a ID disordered chain has a singularity as p'^^^ cx — |e|^^ In^"^ |e| near the 
middle of the band. Notice that the singular increasing of the DOS in a ID system at half-filling is not connected 
with the van Hove singularity, since the latter is located at the band edge for the ID system. Instead it is mediated by 
the interference of impurity scattering and Braggjceflactions at half-filling. Later an analogous singularity has been 
found in the electronic DOS of many ID modelsOLj Since the ID scattering problem is characterized by forward 
and backward scattering processes, ID models with off-diagonal disorder display a singularity at the center of the 
band only if the forward scattering amplitude turns to zero. Small forward scattering with Bragg reflection seems to 
enhance localization and as a result the DOS divergence in the band center is blurred. 

Effects of substitutional impurities on tbe-.nOS of 2D (square, honeycomb and triangular) and 3D (simple cubic) 
lattices have been studied computationallyEZrc3 for cases with diagonal and off-diagonal disorder. The existence of 
a van Hove singularity in the DOS of a two sublattice model for c? > 2 has been shown in [ ^ on the basis of an i 
expansion in disordered systems with n orbitals per site for energies approaching the band center. In all cases diagonal 
disorder has been shown to suppress the van Hove singularity, whereas it is preserved for the 2D simple square lattice 
with off-diagonal disorder. 

We start from the following expression for the DOS 

p('^)(.)=.-|lm|^G«(p,e) (17) 

which expresses p(e) by means of the retarded Green's function G'h(p, e). The new class of diagrams which give the 
dominating contribution to the self-energy of Gji{p, e) is drawn in Fig.^. 
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FIG. 4. Self energy parts which give the leading contributions to the Green's function due to umklapp scattering in a) 
cooperon and b,c) diffuson channels. Diagrams symmetrically conjugate to b) and c) with respect to the impurity line, also 
give a contribution to the self-energy. 



These diagrams represent the contributions to the self-energy in first order in cooperon and diffuson blocks (Ctt (q, e) 
and (k, eV) . Since C-^ (q, e) and Dt^ (k, e) do not carry an external frequency (w = 0) , their expressions are obtained 
from Eqs.([l^) and ( p^ with uj = 0: 
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(18) 
(19) 



As ^ r, both blocks have a diffusion pole for |e|r ^ 1 and ql,kl <C 1. The cooperon and the diffuson blocks 
give rise to logarithmic corrections to the physical parameters for d — 2 and -\/|e|r corrections for d — 3. Therefore, 
diagrams of higher order in the cooperon and diffuson blocks for the self-energy parts are not necessary for 3D systems 
for jejr ^ 1. Nevertheless the logarithmic divergency existing in the first order contribution to the self-energy for 
2D systems requires to examine higher order logarithmic corrections. The structure of the diagrams which contain 
all possible combinations of diffuson and cooperon blocks becomes complicated with increasing number of inserted 
blocks. 

However, as it Jiaspbjeen shown by various theoretical methods, higher logarithmic corrections to the conductance 
cancel each otheitOOH and the leading correction to the conductance is only the first order logarithmic term. 

High order corrections to the self-energy seem to cancel each other also in our case. The second order diagrams 
which give logarithm-squared contributions to the self-energy are drawn in Fig.^. Straightforward calculations show 
that the sum of these diagrams gives zero. 



a) 




b) 
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FIG. 5. Contributions of second order in cooperon and diffuson blocks to the self-energy. The calculations shows a complete 
cancellation between these diagrams. The second-order self-energy diagrams obtained by mutual insertion of the diagrams of 
Fig.^, with a total number of 40, also completely cancel each other at half-filling. 

The retarded Green's function G'fl(e, p) is expressed by the sum of the first-order self-energy parts I](e, p) in Fig.^ 
according to the Dyson equation as 

^ oo 

^-(^'P) - (GO,(,p))-i-E(,p) - E(G?.(e,p))"-(E(.p))" (20) 

By summing the diagrams in Fig.^ one gets the expression for S(e, p): 

E(e,p)= / -^{a(e,k)~— (l--)i?,(e,k)}G?,(-p + k + Qo,e) (21) 

By substituting Eqs.(pO|) and ( |2l| ) into Eq.([l7|) one obtains the following expression for p(e): 

2 



p(6)^p('^)-±Im^A.«S(6) (22) 



n=l 



with 



and 



a,{e)^ir'l ^{^(e, k) - ^(l - ^)D.(6, k)} (23) 



^« = (^/ (0iGr^(.,p)GJ?(-p + po,e) = (-l)"-*-^2.^p('^)!|l^^ (24) 

Here, p^f'' is the DOS for a d-dimensional clean crystal and is given by Eqs.(p^) and ([l^) for d = 2 and d = 3, 
respectively. The sum over n in Eq. (p^) is performed easily using Eq. (p^) : 

^ Vl + a<i(e)(l + yi+M^) ^ + «<i(e) 

where the energy dependence of ad(e) is obtained from Eqs.(p3|), ( p^ and (pj|): 
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(26) 



It can be seen from Eqs.(p5|) and (^6|) that away from half-filling, the effects of Umklapp scattering are weakened 
and ^ Tat « r, as a result of which the quantum corrections to the DOS disappear. In the vicinity of half-filling, 
r « Ttt < Tat and impurity effects become essential. 

The DOS for a 2D system with half-filled energy band can be expressed as 



P^^\^)=p\'^(e) 



(2), 



1 



27reFT7r 



In 



(27) 



4|e|r, 



where p^^^ [t), given by Eq. (p^) , contains a logarithmic singularity in the middle of the band. Eq.(|27|) shows that the 
van Hove singularity in the DOS of a pure 2D square crystal is preserved in the presence of substitutional impurities. 
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However, the central peak becomes narrower due to impurity scattering than that of the van Hove peak in clean 
systems. In the vicinity of the bancL-Genter the energy dependence of p*-^-* (e) is changed from logarithmic dependence 
to the square root of the logarithm,ta 

/'^'H^) = r^(2^^^^'^)i/^lni/2( 1 ) as |e|^0 (28) 
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FIG. 6. Impurity effect on p^^' and p^^' for selected values of the impurity potential. The pure DOS pg^'' and pl^\ corre- 
sponding to T = oo, are also given for comparison. 

In Fig|6^. the dependence of p^^-* (e) on the impurity potential strength is drawn, which displays a narrowing of the 
central peak with increasing disorder strength or decreasing T^r. The DOS of a square lattice with off-diagonal disorder 
computed in [ shows the same effect. Notice that the impurity potential studied in our problem corresponds 

to off-diagonal disorder, since the part of the Hamiltonian Eq.(^ corresponding to the impurity potential can be 
decoupled into diagonal (C/(0)) aad. off-diagonal {U{q) with q 7^ 0) parts. The diagonal part is chosen to be equal to 
zero for the white-noise potential!!^ in the averaging process. 

Unlike 2D systems, substitutional impurities have small effect on the DOS of 3D simple cubic crystals. In the close 
vicinity of half-filling, Eq.(23) for a3(e) is simplified: 
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3 tta/S /n~/^ ■\^ 
— 2"|1 ^ V I'^rll - «j| = "i? + «a/ 
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(29) 



By using Eq.(p9[), on expresses p^^H^) ^ 



(3) 



v/(l + any + aj 
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(30) 



1 



where av{t) and aj{e) are real and imaginary parts of a^{t), respectively. 

Eq.(pO|) shows that a small dip of the DOS is formed at the middle of a half-filled band of a 3D simple cubic lattice 
due to substitutional impurities (see Fig.0b). The depth of this dip increases with disorder strength as 



(3) 



(31) 



The results obtained for the DOS of 2D systems, Eqs.(|25|)-(|2^), and 3D systems, Eqs.(p9|)-(|3l|), show that in contrast 
to ID systems substitutional impurities tend to reduce the DOS on the Fermi surface of 2D and 3D crystals at half 
filling. Indeed, the Dyson singularity, expressing, the enhancement of the states density in the ID lattice with off- 
diagonal disorder, is mediated by impurities ,E3~l3 since the DOS of a regular ID lattice is a smooth function of the 
energy of a half-filled band. 

The mechanism of the relative reduction of the DOS in 2D [Eq.(p7|)] and 3D [Eq.(|30|)] systems is Umklapp scattering, 
the same that increases the ID DOS. However, the nested Fermi surface of a 2D square crystal with nearest-neighbor 
hopping contains also saddle points at {±-^,0} and {0,±^} which cause the van Hove singularity in the DOS. 
Umklapp scattering in this case weakens the central van Hove peak in the DOS, nevertheless it could not damp the 
peak or reverse its sign. We understand the different effects of substitutional impurities on the DOS of ID and higher 
dimensional systems in the following way: The main mechanism of localization in a ID disordered system is backward 
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scattering on impurities .0'EH3 In the process of scattering on an impurity, simultaneous Bragg reflection at half-filling 
reverses backward scattering to forward scattering and vice versa. Therefore for vanishingly small forward scattering 
amplitudeE2lE3 Bragg reflection prevents localization due to impurity scattering. -However, the situation is different 
for 2D and 3D systems. According to the intuitive discussion given by Bergmannj23 localization in 2D systems is due 
to interference of the electronic wave, returning to the starting point after multiple scattering on impurities with only 
a small change of the momentum at each act of scattering, with the wave on the time-reversed path (see also[ ^). In 
this case, the Bragg reflection accompanying the impurity scattering can not destroy the picture of interference and 
as a result can not completely delocalize all states. 

The addition of next-nearest-neighbor hopping terms into the model with strength t' splits saddle points from 
the nested Fermi surface. For the energy spectrum e(k) = —t\cospxa + cospyQ — j cos p^-a cos PyQ + ^] the saddle 
points again lie at {±-^,0} and {0, ±^}. However, the^aptimal nested Fermi surface is realized at t'/t = 0.165 and 
fi — —0.56, with the new nesting vector Qq — 0.91Qo.l3 Our method can be applied also for this dispersion band. 
The singular blocks are again calculated after separating the large momentum transfer Qq and linearizing the energy 

spectrum around the Fermi surface. As a result the same expression (p8|)) for the 2D DOS is obtained. However, pp^' 
in this case has no singularity on the Fermi surface at half-filling. Therefore p'-^-' (e) decreases with the energy around 
the Fermi surface for half filling and vanishes on it. 



IV. CONDUCTIVITY 



The DC-conductivity for the 2D Anderson model at zero temperature has been computed in and the behavior of 
cr as a function of the Fermi energy and the disorder has been studied. However, these numerical results provide limited 
insight into the physical origin of the processes giving contributions to the conductivity. In this section, we study 
low temperature quantum contributions to the conductivity for the model under consideration. Maximally crossed or 
"fan" diagrams are responsible for the quantum interference corrections to the conductivity in the weak localization 
theory.Q Interference between electronic wave functions in the process of multiple scattering on randomly distributed 
impurities changes the mobility of an electronic system and "fan" diagrams give a contribution to the diffusion 
coefficient. The conductivity of simple crystalline systems with substitutional impurities can also be affected through 
changes in the DOS due to Bragg reflection in the scattering processes on the impurities which becomes essential for 
an electronic band close to half-filling. The diagrams which describe first order weak localization corrections to the 
conductivity and represent both the diffusion coefficient and the DOS contributions, are drawn in Fig.0. 



p, e + o; 



p + q,e 




a) 



-p + q + Qo, e + w 
b 




p,e p,e 



c) (c') d) (d') 

FIG. 7. Diagrams giving the first quantum correction to the conductivity, a) and b) characterize the diffusion coefficient 
contributions to cr(tj) from normal and umklapp scattering channels, respectively. E in diagrams c)-d), for the contributions 
from the DOS to o-{uj), is the sum of all self-energy contributions shown in Fig.kl Two other diagrams c') and d') differ from 
c) and d) through the direction of the electron lines. 

Since interference effects are essential for low-dimensional systems, we here will calculate the quantum corrections 
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to the conductivity only for a 2D square lattice with substitutional impurities. The localization problem for a ID 
lattice with off-diagonal impurities has been calculated in [ [l9|-|2^ . 

The quantum correction to the conductivity is calculated according to the Kubo expression 

Ne^ 2e^ f de f (Pp 

aa,f3iuj)^i (5a/3H / — / —-^Vc,{p)vfj{p){G{p,e + uj)G{p,e)) (32) 

muj uj J In J [2tt) 

where G(p, e) is the Green's function and Wq(p) = with a — {x, y) is a component of the electron velocity; the 
bracket (...) denotes averaging over the impurity realizations. 

Far from the half-filled energy band, maximally crossed diagrams with normal scattering are responsible for log- 
arithmic corrections to the conductivity of the 2D weakly disordered non-interacting electron gas.Q These diagrams 
can be redrawn as a ladder series in the particle-particle channel as shown in Fig.^. 

Bragg reflection is intensified as half-filling is approached. In this case, an act of electron scattering on an impurity 
may be accompanied by Bragg reflection on the Brillouin zone boundary. The new diagrams which give a contribution 
to the conductivity as half-filling is approached are drawn in Figs.|^-d'. 

The contribution of the diagrams in Fig. 6a to the conductivity, Saaii^) can be expressed as follows according to the 
Kubo expression Eg. (p^) : 

Saa{ui)^'^J ^ J ^^^CAr(q,e,w)A^(q,e,cj) (33) 

where CAr((q, e,a;) is the cooperon block Eq.(||) for normal scattering and 

j^Vc,ip)vc, i-p + q)G(p, e + u;)Gi-p + q, e + t^)G(p, e)G(-p + q, e) (34) 

For a tight-binding model with nearest-neighbor hopping, the a — {x, y) component of the velocity Wq(p) is given 
as Wq(p) = ta sin(pQ,a). The integration over p can be written as J (2tt)^\v \ f '^^^ then the integral over the energy 
variable ^ is done easily. The van Hove singularity, arising at saddle points of the Fermi surface for the half-filled 



band when |vp| = tay sin Pxa + sin pya vanishes, is removed due to the v^{p) term under the integral over the 

Fermi surface dS. Therefore the value of Aa{q, e, iv) does not strongly differ from that for the free-electron gas model 
and we get 

Ac,{(i,e,Lj) = -pFTHei-e{e + Lj)) (35) 

Using Eqs.(||) and ( |35| ) for GAr(q, e,w) and Aa{q,e,uj) in Eq.(^) one gets 

e'^T / 1 \ 
27r^Tjv ^7 z|w|r|/ 

Far from half-filling and far from other commensurate points the Bragg reflection is weakened and t K±T]y ^ t,^. 
Therefore Saa{(^) in this case represents the conventional weak localization correction to the conductivity.Q Saa{i^) is 
reduced close to half-filling when umklapp scatterings are enhanced and t Tt, tn ■ 

The diagrams which give an essential contribution to the conductivity at half-filling are shown in Figj7|3-@d'. The 
appearance of these new diagrams is due to intensified scattering in the 7r-channel with large momentum transfer. 
The main quantum correction to the diffusion coefficient at half-filling comes from the "fan" diagrams in Fig.|^ with 
the TT-cooperon block. The expression corresponding to this diagram is 

Sabiuj)^— [ -^^-^Crri(l,e,uj)A^{(i,e,uj) (37) 

uj J 2n J {2Try 

where G7r(q, e, uj) is the 7r-cooperon given by Eq.([T^) and 

f d^p 

A^{c[,e,uj) ^ J j^--^Vaip)vai-p + q + Qo)G{p,e + uj)G{p,e) ^^^^ 
G(-p + q + Qo, e)G(-p + q + Qo, e + w) 
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The calculation of Eg. pq) is similar to that for Ai^{q,e,uj). Taking into account the condition that Wa(k + Qo) = 
~VaO^) we obtain the result for ATj{q,e,uj) 

- (1 - 2^rml-2^r\e + ..|) { 1 - + .| ""^^^ + + ^^'^^ + ^^^^ 
Substituting Eqs.Q and @ in Eq.(|7|), the expression for Sab{uj) can be reduced by some simple calculations to 

Ja.H = ^^{(^)^/(^,..r) + ^} (40) 



with 



T . i f\ , I 1 



/( — , Cij = lut) ~ — dx dz 



T^' 2t Jo Jo (z + i)^ (z + uj + i)^ -x + {^) 



2 



^ ^ ' [1 + ^ ]ln(l-icJ + z /l-(— )2 



(41) 



1] ln(l-zcJ-^)} 



The expression for 5ab{^) can be presented for two limiting cases, namely near the middle of the band when r ~ 
T-!T ^ Tjv and far from half-filling when r w r^y <C r^r by using Eq.(^l]) in (^): 

f ln(l - ^ - iujT) + for ^ ^ 1 

^ ' 1 const + i-r^ ( — ) ln2 for ^ w ^ 

So the quantum correction to the conductivity corresponding to the diagram in Figj^b decreases (t(uj) logarithmically 
with the external frequency near half-filling. The contribution of Figj^b vanishes through a small offset from the 
middle of the band. 

Analyzing the effect of substitutional impurities on the conductivity by means of the diffusion coefhcient, which 
is expressed by the diagrams a) and b) in Fig.|^, it can be seen that the logarithmic correction to the localization 
correction to 17(0;) is preserved irrespectively of the band filling. However, the coefficient of the logarithm of the 
quantum correction to the conductivity is changed from in Eq.(p6|) far from half-filling to in Eq.(^) close 
to half-filling. On the other hand, the 7r-scattering mechanism gives a perturbative contribution to the dielectric 
constant e'{i^) ~ lma{uj). The contribution to e'{uj) increases close to half-filling. This also means that there exist a 
few delocalized states at the center of the band. 

The diagrams in Fig.^-d' represent the quantum correction to the conductivity due to changes in the DOS. The 
expression corresponding to the sum of these diagrams can be presented as 



2{5(T,{uj)+5ad{u))^^- j ^ j A«_a(q, e + ^)B(q, 



c.) (43) 



where 



and 



a(q, e) = ^(q, e) - —(1 - -)DM, e) (44) 



B(q, e, 0.) = {(— - 2)e{^e{e + co)) + \ 9(e(6 + to))} (45) 

[1 — 2iT\e + Lu])'' I tn 1 — it\€ + uu\ — iT\e\ ) 

The fact that the contributions from diagrams c' and d' in Fig.^ are equal to those from c and d is taken into account in 
Eq. (p3[) . The cooperon and diffusion insertions into the Green's function are given by Eqs.(p^ and (|l^), respectively. 
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After some routine calculations, Eq.(Q) for the correction to the conductivity due to the DOS is reduced to the form 

2{6a,{u,)+5ad{^)) = -^^3 [l - — (1 - -)] {l{-,^ = ^r) 

^2^^, (l-2^c.r)2 + (^)2 / 

Where I{^,lu — lot) is given by Eg. (p]) . These contributions vanish as expected far from half-filled band, when 
S> T « Tat. However, they give a large contribution to the conductivity near half-filling, r « r^r rjv, giving rise 
to a rapid increasing of the conductivity with the external frequency. So, the total contribution to the conductivity 
is obtained by summing Eqs.(^6|),(|4^), and (^). Far from half-filling, only the normal cooperon gives a contribution 
and, as a result Suaioj) survives: 

5aH = -;^ln(--V) (47) 
ivr^ i\uj\t 

The quantum correction to the conductivity near half-filling is due to the 7r-cooperon and 

= (48) 

with T — > r,r- Although the conductivity decreases with the external frequency, the diffusion coefficient increases with 
approaching half-filling for given w, which means a partial lifting of localization at the center of the band due to a 
few delocalized states. The rapid decrease of the conductivity (Eq.(^)) near half-filling is the result of the impurity 
effect on the DOS. 



V. CONCLUSION 



The character of localization in the 2D disordered electronic gas has been a subject of debate for a long time. The 
question whether delocalized states exist in the center of the band of 2D disordered systems or not is also one of the 
crucial points for the integer quantum Hall effect. 

Searches for delocalized states have been mainly concentrated around the 2D Anderson model. Expectations are 
connected with the existence of the van Hove singularity in the band center of the square lattice, which might givej|iaa 
to a delocalization of sta.ipjpal_half-filling. Various num£ijical approaches have been applied to compute the DOSj23'E2l 
the localization lcngthj3'L3E30'L£l and the conductivityOEa of the 2D Anderson model with diagonal and off-diagonal 
disorder. The computations show that although even a small concentration of diagonal disorder suppresses the van 
Hove singularity, off-diagonal disorder preserves it, however witluajnodified shape. 

The exponential localization of all states for diagonal disordeiH'Ej has been revealed by computational approaches 
whereas -tlie-i existence of quasi-localized states has been predicted for the 2D Anderson model with off-diagonal 
disorder.oE30 In other words, the computation of the effects of off-diagonal disorder provides evidence for the 
transformation from exponential to power-like localization as half-filling is approached. 

In this paper we tried to give a complete picture for the DOS and the conductivity, and also to give some physical 
insight into the processes occuring near half-filling in simple lattices with substitutional impurities. 

The impurity concentration is considered to be small, so that the condition epT ^ 1 or ppl ^ 1 is satisfied and 
the weak localization approximation is applicable. 

Apart from normal scattering with small momentum transfer, umklapp scattering with large momentum transfer is 
shown to be essential near half-filling. This new scattering involves coherent reflection of an electron on the boundary 
of the Brillouin zone for each act of scattering. 

Our diagrammatical approach shows that the new singular blocks, namely the 7r-cooperon C,r(q, Cj^^) and the 
TT-diffuson Z3^(k, e,w), give a large contribution to the DOS and the conductivity at half-filling. The dependence of 
the DOS on the energy and the impurity strength near the middle of the band have been determined by summing 
the leading logarithmically divergent contributions. The results obtained [Eqs.(^7|) and (^)] for the DOS show that 
weak disorder due to substitutional impurities does not remove the van Hove singularity at the center of the 2D band. 
However its energy dependence is strongly changed. The impurity effect on the 3D DOS is a shallow dip on the 
smooth background of the bare DOS in the middle of the band. 
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The effect might be observable in the temperature dependence of the magnetic susceptibihty according to 



X(,)(T) = 2^^% J ^cosh-2(^)p('^)(e) (49) 
Where is the Bohr magneton. Straightforward calculation using Eq.(p7|) for the 2D DOS gives 

X2d{T) (X \n'/-\-^) 

1 

For d = 3, the correction to the magnetic susceptibility can again be calculated according to Eqs.(^9|) and (|3C|). 
However, this correction is negligibly small. 

The first logarithmic corrections to the conductivity of the 2D square lattice with substitutional impurities come 
from both the diffusion coefficient and the DOS. The corrections due to the diffusion coefhcient, Saa{uj) + Sab{uj), 



given by Eqs.(36) and ([12[), decrease in magnitude as half-filling is approached while they remain logarithmically 
dependent on the ext ernal . . fr fi quei icy. Such a partial lifting of localization may be due to a few delocalized states in 
the center of the band.L3E3SL3E3 The imaginary part of cr{uj), which corresponds to the dielectric constant, increases 
close to half-filling. This fact supports the picture of an increase of the electronic mobility at the center of the band. 

Nevertheless, the contribution to cr{uj) due to the DOS suppresses the relative increase in the conductivity coming 
from the diffusion coefficient. 
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